Renormalized Polyakov Loop in the Deconfined Phase of SU(N) Gauge Theory and 
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We use gauge/string duality to analytically evaluate the renormalized Polyakov loop in pure Yang- 
Mills theories. For SU{3), the result is in a quite good agreement with lattice simulations for a broad 
temperature range. 
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INTRODUCTION 

It is well known that a pure SU{N) gauge theory at 
high temperature undergoes a phase transition. This 
phase transition is of special interest because of many its 
aspects can be characterized precisely [1]. In particular, 
the order parameter is given by the Polyakov loop 



L{T) = -tr Pcxp 



^9 



l/T 



dtAn 



(1) 



where the trace is over the fundamental representation, t 
is a periodic variable of period l/T, with T the temper- 
ature, g is a gauge coupling constant, and is a vector 
potential in the time direction. The usual interpretation 
of (1) is as a phase factor associated to the propagation 
of an infinitely heavy test quark in the fundamental rep- 
resentation of the gauge group. 

Until recently, the lattice formulation, still struggling 
with limitations and system errors, and effective field the- 
ories were the main computational tools to deal with non- 
weakly coupled gauge theories. The Polyakov loop was 
also intensively studied (see, for example, [2] and refer- 
ences therein). The situation changed drastically with 
the invention of the AdS/CFT correspondence [3] that 
resumed interest in another tool, string theory. 

In this note we continue a series of recent studies [4-6] 
devoted to a search for an effective string description of 
pure gauge theories. In [4], the model was presented for 
computing the heavy quark and multi-quark potentials 
at zero temperature. Subsequent comparison [7] with 
the available lattice data has made it clear that the model 
should be taken seriously. Later, this model was extended 
to finite temperature. The results obtained for the spa- 
tial string tension [5] and the thermodynamics [6] are 
remarkably consistent with the lattice, too. As is known, 
QCD is a very rich theory supposed to describe the whole 
spectrum of strong interaction phenomena. The question 
naturally arises: How well does the model describe other 
aspects of quenched QCD? Here, we attempt to analyt- 
ically evaluate the Polyakov loop as an important step 
toward answering this question [8]. In addition, a good 
motivation for this test is lattice data revealed recently 



by [9]. 

Before proceeding to the detailed analysis, let us set 
the basic framework. As in [4-6], we take the following 
ansatz for the five-dimensional background geometry 



ds^ = GnmdX^dX"' ^ R^w fdt^ + dx 



:dz' 



(2) 



where = 1/ttT. 5 is a deformation parameter whose 
value can be fixed from the critical temperature [10]. We 
take a constant dilaton and discard other background 
fields. 

In discussing the Wilson and Polyakov loops within 
the gauge/string duality [11], one first chooses a con- 
tour C on a four-manifold which is the boundary of a 
five-dimensional manifold. Next, one has to study fun- 
damental strings on this manifold such that the string 
world-sheet has C as its boundary. In the case of inter- 
est, C is an interval between and l/T on the t-axis. The 
expectation value of the Polyakov loop is schematically 
given by the world-sheet path integral 



(3) 



;L(r))= / DXe 



where X denotes a set of world-sheet fields. is a world- 
sheet action. In principle, the integral (3) can be eval- 
uated approximately in terms of minimal surfaces that 
obey the boundary conditions. The result is written as 
( L{T) ) = '^■ri exp[— S'n], whcrc Sn means a renormal- 
ized minimal area whose weight is Wn- 

CALCULATING THE POLYAKOV LOOP 

Given the background metric, we can attempt to calcu- 
late the expectation value of the Polyakov loop by using 
the Nambu-Goto action for Sw in (3) 



(4) 
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Here Gnm is the background metric (2). In the case 
of interest, this action describes a fundamental string 
stretched between the test quark on C (at z = 0) and 
the horizon at z = z^- Since we are interested in static 
configurations, we choose = i, ^2 = z. This yields 



S 



ttT 



(5) 



where g = A prime stands for a derivative with 

respect to z. 

Now it is easy to find the equation of motion for x 



wfx'/^i + f{x'y 



= 0. 



(6) 



It is obvious that Eq. (6) has a special solution x = const 
that represents a straight string stretched between the 
boundary and the horizon. Since this solution makes the 
dominant contribution, as seen from the integrand in (5), 
we won't dwell on other solutions here. 

Having found the solution, wc can now compute the 
corresponding minimal area. Since the integral (5) is 
divergent at z = due to the factor z~^ in the metric, 
we regularize it by imposing a cutoff e 



5h = 



ttT 



dz w . 



(7) 



Subtracting the i term (quark mass) and letting e = 0, 
we get a renormalized area 



(8) 



•^0 = ;^/ dz{w - — ] + c 



where c is a normalization constant which is scheme- 
dependent. 

Next, we can perform the integral over z. The result 

is 



^o = 0(x/^|Erfi(|)+l-e(^=/^)^)+c. (9) 



In this formula Tc is given by Tc = Vs/"" [5]- 

Combining the weight factor with the normalization 
constant as c = In wq — c, we find 



L(T) = exp 



fll'v^lErfid) 



-1-c 



(10) 



with Erfi(z) the imaginary error function. This is our 
main result. 



NUMERICAL RESULTS AND 
PHENOMENOLOGICAL PROSPECTS 

It is of great interest to compare the temperature de- 
pendence of (10) with other results for the high temper- 
ature phase of SU{N) gauge theory. In doing so, we 
start with lattice QCD. Clearly, = 3 is of primary 
importance. In Fig.l a comparison is shown with the 
recent data of [9]. We see that our model is in a quite 
good agreement with the lattice for a broad temperature 
range 1.05 Tc < 20 Tc. The maximum discrepancy 
occurred at T = 1.05 Tc is of order 15%. It rapidly de- 
creases with temperature reaching 2% at T — 2.2 Tc and 
becoming almost negligible up to 20 Tc. Then, it starts 
to grow back again. 

For completeness, we can fit the value of g to be 0.72 
that significantly improves accuracy. For example, at 
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FIG. 1: The renormalized Polyakov loop in 5'?7(3) gauge 
theory. The solid blue curve corresponds to (10) with g = 0.62 
as fixed from the heavy quark potential at zero T in [7]. The 
dashed green curve represents the "best fit" with g = 0.72. 
In both cases, the value of c is set to 0.10. The dots are from 
lattice simulations of [9]. The red dots are for A''t = 4, while 
the black dots are for A^^ = 8. We do not display any error 
bars because they are quite small, comparable to the size of 
the symbols. 

T = 1.05 Tc it becomes of order 6%. One possible ex- 
planation for the better fit is that we have evaluated (3) 
classically (in terms of strings). If we take into account 
semi-classical corrections, then the value of g gets renor- 
malized. 

For practical purposes, the expression (10) looks some- 
what awkward. Following [6], we expand Sq and L in 
powers of {Tc/TY . If we ignore all higher terms, then a 



3 



final result can be written in two simple forms: 



L{T) « exp 



or 



L(T) «cMl 



T 



(11) 



(12) 



In Fig. 2 we have plotted the results. As can be seen, 
above 2Tc the discrepancy between the expression (10) 
and approximations (11)-(12) is negligible. At lower T 




shown in [15], the deformation parameter 5 of the back- 
ground geometry (2) is tied into the appearance of the 
quadratic corrections. It is not, therefore, surprising that 
we have recovered (11) in our calculations. 

Interestingly, the power law (12) is very similar to that 
observed for the pressure in [16]. Indeed, for T > 1.2 Tc 
the pressure is simply p/T* « /p„.(l - (Tc/Tf). 



CONCLUSIONS 

In this note we have evaluated the Polyakov loop us- 
ing the now standard ideas motivated by gauge/string 
duality. A key point is the use of the background metric 
(2) which is singled out by the earlier works [4-6]. (Note 
that there is no need for any free parameters except a 
scheme-dependent normalization constant c.) The over- 
all conclusion is that the same background metric results 
in a very satisfactory description of the Polyakov loop as 
well. Of course, we still have a lot more to learn before 
answering the question posed at the beginning of this 
note. 
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FIG. 2: A comparison of different L{T) curves for SC/(3) 
gauge theory. As in Fig.l, the solid blue curve corresponds to 
(10) and the dots are from lattice simulations of [9]. The blue 
dashed curve corresponds to the exponential law (11). The 
black dot-dashed curve corresponds to the power law (12). In 
all the cases, g = 0.62 and c = 0.10. We display error bars 
only if they are comparable to the size of the symbols. 

the approximation (11) (exponential law) is poor. It 
shows a significant deviation from the lattice. In par- 
ticular, the discrepancy occurred at T = 1.05 Tc is of or- 
der 27%. On the other hand, the agreement between the 
approximation (12) (power law) and the lattice is spec- 
tacular. For the temperature range 1.05 Tc < T < 20 Tc 
the power law provides a reliable approximation to lat- 
tice QCD with accuracy better than 5%! Moreover, one 
can use it to describe all available lattice data of [9] at 
lower T. Then, the maximum discrepancy occurred at 
the lowest available value T = 1.012 Tc is of order 7%. 

It is worth noting that the exponential law has been 
suggested in [12] based on a dimension- two condensate 
{A^) [13]. Such a condensate as well as its possible links 
to the UV renormalon and l/Q^ corrections got inten- 
sively discussed in the QCD literature [14]. As was first 
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